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We present a microscopic calculation of the complete quasiparticle interaction, including cen- 
tral as well as noncentral components, in neutron matter from high-precision two- and three-body 
forces derived within the framework of chiral effective field theory. The contributions from two- 
nucleon forces are computed in many-body perturbation theory to first and second order (without 
any simplifying approximations). In addition we include the leading-order one-loop diagrams from 
the N 2 LO chiral three-nucleon force, which contribute to all Fermi liquid parameters except those 
associated with the center-of-mass tensor interaction. The relative-momentum dependence of the 
quasiparticle interaction is expanded in Legendre polynomials up to L — 2. Second-order Pauli 
blocking and medium polarization effects act coherently in specific channels, namely for the Landau 
parameters /i, ho and go, which results in a dramatic increase in the quasiparticle effective mass 
as well as a decrease in both the effective tensor force and the neutron matter spin susceptibility. 
For densities greater than about half nuclear matter saturation density po, the contributions to the 
Fermi liquid parameters from the leading-order chiral three-nucleon force scale in all cases approxi- 
mately linearly with the nucleon density. The largest effect of the three-nucleon force is to generate 
a strongly repulsive effective interaction in the isotropic spin-independent channel. We show that 
the leading-order chiral three-nucleon force leads to an increase in the spin susceptibility of neutron 
matter, but we observe no evidence for a ferromagnetic spin instability in the vicinity of the satu- 
ration density po- This work sets the foundation for future studies of neutron matter response to 
weak and electromagnetic probes with applications to neutron star structure and evolution. 
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I. INTRODUCTION 

In a series of recent articles [l|, Q , we have revisited the Fermi liquid description of infinite nuclear matter in 
the context of modern two- and thrcc-nuclcon interactions derived within the framework of chiral effective field 
theory. In these studies we limited our discussion to the central components of the quasiparticle interaction in 
a medium of spin-saturated symmetric nuclear matter characterized by the nucleon density p — 2k^/3ir 2 . In 
the vicinity of the saturation density po, the central components of the quasiparticle interaction are strongly 
constrained by the properties of bulk nuclear matter and its low-energy excitations. In refs. 0, it was found 
that microscopic calculations of the quasiparticle interaction within many-body perturbation theory yield an 
accurate description of the nuclear matter compressibility, isospin asymmetry energy and spin-isospin response 
only with the inclusion of leading-order medium effects, which arise at second order in a perturbative calculation 
with two-body forces and at first order in a calculation with three-nucleon forces. 

In the present paper we extend these calculations to pure neutron matter, which we treat as a normal 
(non-supcrfluid) Fermi liquid (the generalization of Fermi liquid theory to superfluid Fermi systems has been 
performed in refs. 041] )■ Such a quasiparticle description has been used in previous works to understand 
neutrino emissivity in neutron stars [6|-|9( as well as the spin response of neutron star matter to strong magnetic 
fields [Tol - [l3j . In such calculations, the Landau parameters that characterize the quasiparticle interaction 
have been computed with microscopic two-body nuclear forces or with phenomenological Skyrme and Gogny 
effective interactions. Qualitative differences arise between the predictions of microscopic and phenomenological 
forces, particularly with respect to the possibility of bulk magnetization or even the existence of a spontaneous 
ferromagnetic phase transition of neutron matter at several times nuclear matter saturation density fl3l - [la) . 
As discussed in Section [Vl three-neutron forces give rise to Fermi liquid parameters that scale approximately 
linearly with the neutron density (at leading order), and one of the primary goals of the present work is to 
better understand the role of three-neutron correlations in describing the dynamical response of neutron matter 
to weak or electromagnetic probes. 

Both the magnetic susceptibility of dense neutron matter as well as neutrino elastic scattering, absorption and 
emission rates are sensitive to the inclusion of noncentral components of the quasiparticle interaction 8, 9|, ITl| , 
which have been introduced in refs. [la. Il7j. Such interactions include (in the long- wavelength approximation) 
the exchange tensor interaction, proportional to Si2(q) = 3<7i • qa2 ■ q — a% • (?2, which appears already in 
the free-space nucleon-nucleon potential where it is dominated by one-pion exchange, as well as center-of-mass 
tensor and cross- vector interactions, proportional to Si2(P) and (c?i x a^) • (q x P) respectively. In these spin- 
dependent operators, q — p\ — P2 is the relative momentum and P = pi + P2 the center-of-mass momentum of 
the two quasiparticles on the Fermi surface ( |_Pi.2 1 = &/)• When including effects of the medium in the form 
of loop integrals over the filled Fermi sea of neutrons, all noncentral interactions can be generated. In fact, 
the second-order contributions from two-body forces produce exchange tensor, center-of-mass tensor and cross- 
vector interactions. Note that the spin-nonconserving cross-vector interaction can only be generated through 
polarization (particle- hole) corrections to the effective interaction and not by ladders [TtJ ■ To date there has been 
no study of three-nucleon force contributions to either the center-of-mass tensor or cross-vector interactions, 
while in ref. [l8T | loop diagrams involving intermediate A-isobar excitations were analyzed and shown to generate 
an exchange tensor interaction in symmetric nuclear matter. 

The present work sets the foundation for a systematic study of the role played by second-order corrections and 
three-nucleon forces in determining the noncentral components of the quasiparticle interaction. Their effect on 
neutron star structure and evolution will be studied in future work. We employ many-body perturbation theory 
and compute all contributions to the quasiparticle interaction for the case of two neutrons on the Fermi surface 
interacting in a background medium of pure neutron matter. We employ high-precision chiral two- and three- 
nucleon interactions whose unresolved short-distance components are encoded in a set of contact couplings 
proportional to low-energy constants fit to elastic nucleon-nucleon scattering phase shifts and properties of 
light nuclei [l9l - [2l| . To improve the convergence of the microscopic calculation of the quasiparticle interaction 
in many-body perturbation theory, we employ renormalization-group methods for integrating out momenta 
beyond the scale of A ~ 2.0 fm -1 , which results in the nearly universal two-nucleon potential Mow-k [H, |23j . 
The contributions to the quasiparticle interaction from the leading-order chiral three-nucleon force (with scale- 
dependent low-energy constants) is computed to one-loop order. 

The present paper is organized as follows. In Section [TT] we give a brief review of Landau's theory of normal 
Fermi liquids and discuss the microscopic approach for computing the quasiparticle interaction. We then present 
a general method that can be applied to any free-space nucleon-nucleon force given in a partial-wave basis for 
extracting the scalar functions that multiply the various spin-dependent operators (namely the central, spin- 
spin, exchange tensor, center-of-mass tensor and cross-vector terms) occuring in the quasiparticle interaction. 
We then discuss the microscopic origin of the spin-nonconserving cross-vector interaction, which at second- 
order results exclusively from the interference of a two-body spin-orbit force with (in principle) any non-spin- 
orbit component in the free-space nucleon-nucleon interaction. In Section IIIII we test the intricate (numerical) 
calculations of the second-order contributions to the quasiparticle interaction by means of model interactions 
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that can be solved partially analytically. In Section lPVl we present analytical formulas for the Landau parameters 
of the quasiparticle interaction arising from the leading-order N 2 LO chiral three-neutron force. The results of 
the corresponding calculations with two- and three-neutron forces at several resolution scales are presented and 
discussed in Section [V] We end with a summary and conclusions. 

II. QUASIPARTICLE INTERACTION IN NEUTRON MATTER 

A. General structure of the quasiparticle interaction and spin-space decomposition 

In Landau's theory of normal Fermi liquids [24l. [25j. the quasiparticle interaction J r (piS\ti;p2S2t2) derives 
from the change in the total energy density due to second-order variations in the particle occupation densities 

d£ = y^epfapst + - ^ J 7 {piSiti;p2S 2 t2)Snp lSl t 1 Snf 2S2t2 H , (1) 

pst pisiii 

P2S 2 *2 

where Si and ti label the spin and isospin quantum numbers of the zth quasiparticle. The most general form 
for the effective interaction between two quasiparticles in pure neutron matter in the long-wavelength limit is 
given by |17| 

F{pi,P2 ) = f(pi,P2 ) + g{p\,p2 )&\ ■ VI + h(Pl,P2 )Sl2(q) + k(pi,p 2 )Sl 2 (P) 

+l(pi,p 2 )(a 1 Xa 2 )-(qxP), (2) 

where q=pi —p 2 is the momentum transfer in the exchange channel, P = p\ +P2 is the conserved center-of-mass 
momentum and the tensor operator Si 2 (v) is defined by Si 2 (v) — 3cti ■ va 2 ■ v — 5\ ■ B 2 . The interaction in eq. 
(|2|) is invariant under parity and time-reversal transformations as well as under the interchange of the particle 
labels. However, due to the presence of the medium, Galilean invariance is no longer manifest, leading to new 
operator structures (namely Si 2 (P) and A\ 2 {q, P) — (ai x a 2 ) ■ (q x P)) that depend explicitly on the center- 
of-mass momentum P. Since the two quasiparticle momenta lie on the Fermi surface \p±\ — \p 2 \ — kf, the 
remaining angular dependence of the quasiparticle interaction is conveniently expanded in Legendre polynomials 
of cos 6 =pi ■ p2 ■ 

oo 

X(P1,P2) = J2 XL(k f )P L (co S 9), (3) 

where x represents /, g, h, k, or I, and the angle 9 is related to the relative momentum q = \pi~p 2 \ through the 
relation q = 2kf sin (0/2). 

The Fermi liquid parameters of nuclear matter can be either extracted from experiment (due to direct relations 
between particular Landau parameters and nuclear observables [25l |26| ) or computed microscopically within 
many-body perturbation theory (27| . Given a strongly interacting normal Fermi system at low temperatures, 
it is not clear that a pertubative approach is justified. However, in the case of nuclear and neutron matter, 
there are strong indications that renormalization group techniques [22t l23l |28| may render nuclear interactions 
perturbative for a wide range of densities when evolved down to a cutoff scale A < 2.1 fm -1 . In the framework 
of many-body perturbation theory, the quasiparticle interaction is extracted by functionally differentiating 
the contributions to the ground-state energy density twice with respect to the particle occupation numbers. 
Previous work [l|, [2] has shown that a satisfactory description of bulk nuclear matter properties around the 
saturation density can be obtained with chiral and low-momentum nuclear interactions only with the inclusion 
of the leading-order Pauli-blocking and medium-polarization effects from two- and three-nucleon interactions. 
This requires a second-order perturbative calculation in the case of two-nucleon interactions together with the 
first-order perturbative contribution from three-nucleon forces. These systematic calculations suggest that a 
consistent microscopic description of the quasiparticle interaction in neutron matter, for which there is much 
less empirical information, can be achieved. 

The relations between specific Fermi liquid parameters and bulk properties of the medium are well known 
[25|,[26|. The compression modulus of neutron matter (at density p = k^/3iT 2 ) 

K=^(l + *b), (4) 

is related to the Fermi liquid parameter F = N f (with iVo = kfM* /tt 2 the density of states at the Fermi 
surface) which represents the isotropic part of the spin-independent quasiparticle interaction. Note that the 
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compressibility JC — k'j d 2 E / dkj + 4kf dE jdkf is determined by both the curvature and slope of the energy per 
particle E. In eq. ((4]) the quasiparticle effective mass M* is given by 



M* _ F 1 



(5) 



with M n = 939.6 MeV the free neutron mass and F\ — Naf\. Considering only the central components of the 
quasiparticle interaction, the neutron matter spin susceptibility is given by 



X = AV, 



N 



1 + Go 



(6) 



where fi n = —1.913 is the free-space neutron magnetic moment (in units of the nuclear magneton) and Go = 
Nogo. The presence of noncentral components in the quasiparticle interaction that couple quasiparticle spins 
to their momenta results in effective charges (magnetic moments) that are not scalars under rotations of the 
quasiparticle momentum. The resulting expression for the spin susceptibility then involves both longitudinal 
and transverse components of the magnetic moment [T3, [H| • 

The first calculation to extract all components of the quasiparticle interaction given in eq. ^ for an arbitrary 
two- neutron force was performed in ref. [17j . In the following we present a general method to project out the 
various momentum-dependent scalar functions /, g, h, k and / of the quasiparticle interaction in eq. @. This is 
achieved by taking specific linear combinations of the spin-space matrix elements of the quasiparticle interaction. 
The form of this matrix will of course depend on the choice of coordinate system. For q = p\ — p% = qe z and 
P = P1+P2 = P&x, the spin-space matrix elements (m s \F(pi, p2)\m' s ) in terms of the scalar functions /, g, h, k, 
and I are given by 
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The upper left 3x3 submatrix gives the nine triplet-triplet matrix elements, while the fourth row and fourth 
column give the matrix elements that include the singlet state. In this coordinate frame, the exchange tensor 
interaction h(pi,p2 )Si2(q) gives nonzero contributions only for the diagonal triplet spin-space matrix elements. 
The center-of-mass tensor force fc(f?L,2?2 )S±2(P) contributes to the triplet diagonal matrix elements as well 
as the matrix elements mixing m s — ±1 with m' s = =pl. The cross-vector interaction l(pi,p2 )^4i2(<7, P ) is 
nonvanishing only in spin-nonconserving transitions with \m s — m' s \ = 1. We note that it is not possible to 
separate the exchange tensor from the center-of-mass tensor interaction when considering only diagonal spin- 
space matrix elements, since they appear in the same linear combination 2h — k. An alternative choice of 
coordinate system in which q = qe x and P = Pe z just interchanges the contributions of the exchange tensor 
and center-of-mass tensor operators (and I comes with opposite sign). 

From the above matrix, the scalar functions multiplying the spin-dependent operators are extracted as the 
following linear combinations of the spin-space matrix elements: 

J=(2^ 1+ ^ +J- s i0 )/4, 
g =(2Tt A +JS i0 -37S i0 )/12, 
h = (Fl 1 +Fl_ 1 -Fl )/6, 
k=(F t h _ 1 )/3, 

l = (Fi%)/V2, (7) 

where the superscripts s,t distinguish spin-singlet and spin-triplet states of two quasiparticles, and the two 
subscripts label m s and m' s , respectively. Only the function l(pi,P2) multiplying the cross- vector operator 
Ai2(q,P) depends on a singlet-triplet mixing matrix element. 



B. First-order contribution 



With the relations given in eq. ([7]) we need to compute only the quasiparticle interaction in different total 
spin states. The first-order perturbative contribution to the energy density is given by 



S2N = n y^(fci s i; k 2 S2\V\kisi; k 2 s 2 )n 1 n2 7 



(8) 
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FIG. 1. Diagrams contributing to the second-order quasiparticle interaction (all interactions represented by wavy lines 
are antisymmetrized): (a) particle-particle diagram, (b) hole-hole diagram, and (c) + (d) particle-hole diagrams. 

where V = (1 — P\ 2 )V denotes the antisymmetrized two-body potential, rij = 9(kf — \kj\) is the usual zero- 
temperature Fermi distribution and the summation includes both spin and momenta. The corresponding con- 
tribution to the quasiparticle interaction reads 

T^>(p lSl ;p 2 s 2 ) = (piSi;p 2 s 2 \V\piSi;p2S 2 ) = (12|F|12). (9) 

Setting the relative momentum q — p\ — p 2 along the e 2 -direction and projecting onto Legendre polynomials 
Pl{Pi 'P 2 ), the Fermi liquid parameters are obtained from eq. (J9j> in terms of the matrix elements of V depending 
on p = q/2 in a partial-wave representation: 

F L {k f ; Sm s m' s ) = 2(2L + 1)^2 ^ ( X + M)' +5 ) V( 21 + l )( 21 ' + 1) Sm,\ JM) 

IV J 

x{l'0Sm' s \JM) f ' dpj^{plSJM\V\pl'SJM)P L {l - 2p 2 /k 2 f ), (10) 
Jo / 

where we are following the normalization convention in ref. The leading-order contribution is simply a 
kinematically-restricted form of the free-space interaction, which contains no center-of-mass dependent com- 
ponents. Therefore, the only nonzero contributions to the quasiparticle interaction are /, g, and h. This is 
further reflected in the partial wave decomposition through the obvious property M = m s = m' s implied by the 
Clebsch-Gordan coefficients in eq. (|TU)) . 



C. Second-order contributions 



The contribution to the energy density from the two-neutron force at second-order in many-body perturbation 
theory has the form 

5(2) _ 1 ^ |(12|y|34)| 2 Ti in2 (l-n3)(l-n4) 

where in a plane- wave basis ej = k^/2M^*' is the single- particle energy associated with the momentum kj. 
Functionally differentiating twice with respect to the quasiparticle distribution functions yields four different 
contributions to the quasiparticle interaction: 

|(12|t/|34)| 2 (l-n 3 )(l-«4) , l r |(12|F|34)| 2 n 3 n 4 



T (2),^ , - . v ly 12T/ 34 )Y l-»3 l-» 4 , 1^ 



2 ^ ei + £2 - £3 - £4 2 ^ £ 3 + £ 4 - £l - £ 2 

y |(13|V"|24)| 2 n 3 (l - n A ) y |(13|V"|24)| 2 n 4 (l - n 3 ) 
~^ £l + £3 - £2 - £4 4^ £1 + £4 - £2 - £3 

A graphical representation of these terms is given in Fig. [T] The first two terms in eq. (fT2|) are called the 
particle-particle (pp) and hole-hole (hh) contributions, which have a very similar structure. The particle-particle 
contribution is given by 

ly f d 3 k 3 d 3 fc 4 (pip2Sm s \V\k 3 k 4: Sm s )(k 3 k4 : Sm s \V\piP2S'rn' s ) 



J™(p 1 p 2 -,Sm 3 S'm' s ) = ~J2 J 



(2tt) 3 (2tt) 3 £ pi + £ P2 - £fc 3 - £fe 4 

OTfbg 

x (1 - n 3 )(l - n i ){2^)H{p l +p 2 -k 3 - fc 4 ), (13) 
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where we allow for the possibility that m s ^ m' s . The hole- hole contribution is easily obtained from eq. (fT5|) by 
changing the sign of the energy denominator and replacing the two particle distribution functions (1 — 713) (1 — 714) 
with hole distribution functions 77,3774. In computing the pp and hh diagrams, we found it more convenient to 
align the total momentum P along the e^-direction and the relative momentum q along the e^-direction. For this 
choice of coordinate system h and k are just interchanged in eq. ([7]). In a partial- wave basis the particle-particle 
contribution in eq. (|13[) reads 

or 1 1 r k F r<x> 

J** (Sm s m' s ) = V dpp dkk 2 N{l l ml 2 inh'm'l i fh)P^{Q)P^ (0) 

4?r k F , 7o J P 

fhrh s J J M 

M , , /•min{a:o,l} 

v JU " Ah+h-h-UnJM r JM r J'M r J'M / > p m i \ pin / \ 

x o 1.2 L -I 1 mSm J t l2mSm I t I 3 m'Smi t '!4iaSm, / ax r l 2 \ X ) r li \ X ) 

P K Jmax{-ij,-l} 

x{phSJM\V\kl 2 SJM){kl4,SJ'M\V\pl 3 SJ'M)P L (l - 2p 2 /k 2 F ), (14) 

where P, m are the associated Legendre functions, x — cosdk, p = (pi — P2)/2, k — (k 3 — £t)/2, xq = (fc 2 — 
p 2 )/(2ky/k 2 F -p 2 ) and A? \hml 2 rhl z m' 7 4 m) = N™N™N™' N™ with AT/" = ,/(2l + 1)(Z - m)!/(Z + m)!. For 
777 s = 777,' s , this expression agrees with eq. (21) in ref. From the underlying parity invariance of the two- 
neutron interaction, which preserves the total spin 5, the pp and hh contributions cannot give rise to the 
spin-nonconserving cross-vector interaction. However, spin-orbit and tensor terms in the free-space neutron- 
neutron interaction do not conserve m s and therefore can, at second-order in the pp and hh diagrams, give rise 
to a center-of-mass tensor interaction. This mechanism is exemplified in Section IlIII with spin-orbit and tensor 
two-body model interactions. 

We split the second-order particle-hole (ph) contribution into two pieces J 7 ^ + (see Fig. [1} . The first 
term arises from a coupling of the incoming or outgoing quasiparticle 1 to a hole state, while for the second 
term quasiparticle 1 couples to a particle state. The explicit expression reads 

( c ) _ _ , f d 3 h d 3 k 4 (pik 3 S 1 S 3 \V\p2k4S2S4}(p 2 k4s' 2 S4\V\p 1 k 3 S l 1 S 3 ) 



(2tt) 3 (2tt) 3 e P2 + e ki - e Pl - t H 

x 77 3 (1 - 77 4 )(27r) 3 (5(pi + k 3 -p 2 - fc 4 ), (15) 

where we allow for the possibility that the spin states of an incoming and outgoing quasiparticle can be different. 
In contrast to the treatment of the ph contribution to the central components of the quasiparticle interaction 
described in ref. [l|, here we have found it convenient to set p\ — p2 = q e z and p\ + p 2 = Pe x . In this case we 
write 

9 T -X- 1 p2tz pkf pi 

J r ^ (c) (sis 2 si4) = 3 / dcos0P L (cos0) d(j) 3 dk 3 k 2 3 dcos8 3 



(p'hm 1 sis 3 \V\k'l2m2S2S4){k'l4m4s'2Si\V\p l l 3 Tn 3 s' l s 3 ) 



x E 

x cos((m 3 - mi + 777 2 - bk)^)^ 1 (cos 6 p ,)P™ 2 (cos 6 k ,)P™ 3 (cos <V )C 4 ( cos d v ) 

xi l2+l3 - h - h N(l imi l 2 m 2 hm 3 hm 4 )— , ^" 2 , (16) 

kfk 3 cos 3 sin 6/2 + k j sin 0/2 

where p' = [p 1 -k 3 )/2, k' = (p 2 ~k 4 )/2, y = fc/(l-2sin0/2) and z = (fc 2 - fef - 4fc 2 sin 2 0/2)/(4fc/fc 3 sin0/2). 
The product of exponentials coming from the spherical harmonics has been simplified to a cosine by noting 
that the imaginary part vanishes and that <fi p > = cf>k' ■ The expression in eq. (|16j) can be further written out in 
terms of partial-wave matrix elements by first coupling to total spin S and then total angular momentum J. 
The extraction of the scalar functions /, g, h, k and I is achieved by taking appropriate linear combinations of 
the sixteen spin-space matrix elements J r L(sis 2 s' 1 s' 2 ), namely, 

(Sm B \jf\S f m f a )= J2 cS CX> C t™l sjt^A)- (17) 

L / 2 x 2 2 2 1 2 2 

SiS 1 S2S 2 

It is a good check of the calculation that the resulting spin-space matrix on the left-hand side of eq. (IT71) is of 
the form introduced in subsection III Al 

T he p article-hole polarization contribution can give rise to all noncentral interactions, as pointed out first in 
ref. [17|. However, the microscopic origin of the spin-nonconserving cross- vector inter ation should be clarified. 
In fact, neither tensor forces nor spin-orbit forces alone, when iterated in the particle-hole channel, generate the 



7 



cross- vector interaction. This can be seen from the spin structure of eq. (|15l) . In order to produce a nonvanishing 
singlet-triplet mixing matrix element, we can consider without loss of generality the spin- flip transition 

(n n) ~ J2^ 3 t s 3\y\Ah t s 4 )(p 2 k 4 1 S4 imfc 3 1 s 3 ). (is) 

S3S4 

It is easily shown with momentum conservation (pi + k 3 = P2 + fei) that p' — k' = p\ — §1 = qe z . In this case, 
the free-space tensor force (~ Serf erf — &i • 02) vanishes for |Am a | = 1, and therefore one of the two matrix 
elements in eq. (|18[) will vanish for any values of S3 and S4. Similarly, for a free-space spin-orbit interaction of 
the form iV so (<7i +ct 2 ) • (p' x k' ), the vector p' x kl lies in the x — y plane and in this case the spin-orbit operator 
(c?i + <72 ) • (p 1 x k! ) has very restricted matrix elements which are nonvanishing only in spin-triplet states with 
I Am^ I = 1- Again, one of the matrix elements in eq. (|18[) will vanish for any possible values of S3 and S4. These 
arguments indicate that at second-order only the interference of a spin-orbit interaction with any other (non- 
spin-orbit) component can produce the cross- vector interaction. In the following section, we will demonstrate 
that in fact central, spin-spin and tensor components all give nonvanishing interference terms. Finally, we point 
out that although individually J 7 ^ and can have forbidden matrix elements for |Am s | = 1 in the triplet 

states, when summed together such terms cancel exactly. All other spin-space matrix elements of jty and J^j- 
are identical. These features serve as a good check of the involved calculation of the quasiparticle interaction 
at second order. 



D. Contribution from chiral three- neutron forces 

Finally we consider the leading-order contribution from the N 2 LO chiral three-nucleon force [2(| in pure 
neutron matter. Only the components of the two-pion exchange three-nucleon force proportional to the low- 
energy constants c\ and C3 remain for neutrons: 

with parameters gA = 1.29, /„■ = 92.4 MeV and m w — 138 MeV (average pion mass). The quantity <fj is the 
difference between the final and initial momenta of neutron i. In the following we employ two different choices 
for the low energy constants c\ and C3 in eq. (|T9| . When combined with the bare chiral N 3 LO nucleon-nucleon 
potential [l9| we choose the values ci = —0.81 GeV -1 and C3 = —3.2 GeV^ 1 , whereas with the low-momentum 
nucleon-nucleon potential Vi ow -k we use ci = — 0.76GeV _1 and C3 = — 4.78GeV _1 0, [3(j. This variation in 
the low-energy constants (as well as the resolution scale A) provides a means for assessing theoretical uncertainty 
at a given order in many-body perturbation theory. 

The first-order contribution to the energy density of neutron matter has the form 

1 r V f d^h^kjd^ki - 

where V% n = V% n (l — P12 — P23 — P13 + P12P23 + P13P23) is the fully antisymmetrized three-neutron interaction 
and rij = 6(kf — \kj\) + (2ir) 3 5 3 (kj — p*j)5np. a .. Functionally differentiating twice with respect to the two 
quasiparticle distribution functions yields 

^(PUP2) = J ^ ni {il2\V 3n \il2). (21) 

As we will see in Section ITVl the form of the N 2 LO chiral three-neutron force is sufficiently simple that in most 
cases analytical formulas for the Landau parameters are possible. We note that this leading-order contribution 
to the quasiparticle interaction from the N 2 LO chiral three-neutron force is nearly equivalent to the effective 
interaction calculated previously in refs. [3lTj33| . although the quasiparticle interaction represents a restricted 
kinematical configuration for which the two interacting particles lie on the Fermi surface \pi^\ = kf. 



III. BENCHMARK CALCULATIONS WITH MODEL INTERACTIONS 



In order to verify the spin-decomposition techniques and the accuracy of the intricate numerical calculations 
involved in the second-order calculation of the quasiparticle interaction, it is useful to examine simple model 
interactions that are amenable to (partial) analytical solutions. Intermediate-state momentum integrations and 
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FIG. 2. Diagrammatic contributions to the second-order quasiparticle interaction in neutron matter. Diagrams (a)-(d) 
correspond to the sum of particle-particle and hole-hole contributions, while (e)-(g) together comprise the particle-hole 
contribution. Medium insertions are denoted by the short double lines, and the labels (b) and (d) refer to the crossed 
terms of (a) and (c). Reflected diagrams are not shown. 



spin traces are carried out explicitly without decomposing the interaction into partial waves. For details on this 
approach to computing Fermi liquid parameters, see ref. [18] . The diagrammatic contributions are shown in 
Fig. [21 where the double dash on a fermion line represents a "medium insertion" . It is defined as the difference 
between the free-space propagator and the in-medium propagator: 

gOgl - kf) ejkf - \p\) 



p -p 2 /2M n + ie p -p 2 /2M n -ie / 

= -2 /oaf T ■ - 2n5( Po -p 2 /2M n )e(k f - \p\). (22) 

po -p 2 /2M n + le 

Comparing with the expression given in eq. (fT2"j) . the sum of diagrams (a)-(d) corresponds to Fp P + IF). J. 
Expanding (1 — n^){\ — 714) — 1 — 713 — m + 713774 in the pp diagram, we see that the term proportional to 
713714 cancels the hole-hole contribution. Hence, the sum of the particle-particle and hole-hole diagrams gives 
just a free-space contribution and two terms with one medium insertion. The remaining diagrams (e)-(g) in 
Fig. [5] correspond to the particle- hole contribution in eq. (|12[) . In the following we consider only the Fermi 
liquid parameters hf,(kf), kf,(kf), and lh(kf) associated with the noncentral components of the quasiparticle 
interaction. Detailed calculations for the central components have been presented in ref. [l[ . We provide explicit 
formulas for the relevant Landau parameters up to L = 1 obtained by first decomposing the effective interaction 
into the relevant operators and then projecting the expansion coefficients onto Legendre polynomials Pl(pi ■P2)- 



A. Pseudoscalar boson exchange to second order 

We begin by considering pseudoscalar boson exchange with a "form factor" modification: 

g 2 

V(q) = - -, 2 ai-qa 2 -q, (23) 

where g is a dimensionless coupling constant and m is the mass parameter chosen to be sufficiently large to 
achieve good convergence in momentum integrals and partial wave summations. The momentum transfer is 
denoted by q, and the squared denominator in eq. (1231) insures that all loop integrals converge. As discussed 
previously, the leading-order (free-space) contribution to the quasiparticle interaction contains only the exchange 
tensor interaction (in addition to the central components). Its L — 0, 1 projections are given by 

^ (1) = ^{i^ 1 + 4u2 )-T^}' ^ 



= &A l 4^- Hi + 4. 2 ) - i + -A-i \ . (25) 



7 2 

i 2 \ 4it 4 ' u 2 l + 4u 2 

with the dimensionless parameter u = kf/m. At second order, the direct contributions J 7 ^) and T^ 2 ^ have no 
noncentral components. The crossed term (b) of the iterated pseudoscalar-exchange diagram gives a contribution 
to the exchange tensor interaction: 

, (t s(2f,) 9 4M n f 1 , l + 2u 2 1 



2 



487rm 3 [2u 2 1 + u 2 l + 2u 

/•", 1+ 4a; 2 +8x 4 , , 

+ / ax — — (arctan2x — arctanz) > , (26) 

J Q u 2 (1 + 2x 1 y 
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Modified pseudoscalar boson exchange (fc/ = 1.7fm 1 ) 




h [fin*] 


feo [fin 31 ] 


hi [fm*] 


fci [fm 2 ] 


2nd(pp) 


0.134 


-0.142 


-0.166 


0.053 


2nd(hh) 


0.027 


-0.061 


-0.057 


0.074 


2nd(ph) 


1.394 


-0.001 


0.912 


-0.001 


Total 


1.555 


-0.205 


0.689 


0.125 


Analytical 


1.552 


-0.209 


0.719 


0.126 



TABLE I. The L = 0, 1 noncentral Fermi liquid parameters from the pseudoscalar exchange interaction in eq. (I23p at 
second order. We compare the sum of the particle-particle, hole-hole and particle-hole diagrams computed numerically 
to the semi-analytical results of eqs. (|26[) - (|37[) . 



hi(k f ) 



(26) 



g 4 M n f 



167TTO 3 |_ 1 + 2u 2 



2u 4 



In- 



1 + u 2 



u u 2 - 2x 2 

+ I dx—r, ^-(1 + 4a; 2 + 8a; 4 )(arctan2a; — arctana:) 

u 4 {\ + 2x 2 ) i Jy ' 



(27) 



Iterated pseudoscalar exchange does not include medium modifications and therefore cannot generate an effective 
interaction that depends explicitly on the center-of-mass momentum P. In contrast, the crossed term (d) from 
the planar box diagram with Pauli blocking gives rise to an exchange tensor force and center-of-mass tensor 
force (as noted in Section HH the particle-particle and hole-hole diagrams cannot generate the cross-vector 
interaction). The associated L = 0, 1 exchange tensor Fermi liquid parameters are 



h {kf) 



(2d) 



9 4 M n 
I2ir 2 k) Ju 



i ,i ,i 
dx dy I dz ■ 



-l J-i 



A] 2 [u- 2 +B] S 



x < yz — x(x + y + z) — 



4x 2 \y + z\ 
A + B-4 



+ (1 _ x2)2 (A + B) v 0(W) 



(A + B - 4)VW 



(28) 



4 2 

x — x yz 



9 4 M n [\ f 1 f\ x 2 
^ 2 kj J Q dX JJ y ]J Z [u- 2 +A] 2 [u- 2 +B] 2 

+x(y + z)(2 yz) + I(3»V - y 2 - z? 1) + 

2^2, s (A + B) V 6(W) 

+x(l - arr(y + z- x)- ' K ' 

while the center-of-mass tensor Fermi liquid parameters are given by 



k (k f ) 



(2d) 



9 4 M n 
12ir 2 k 3 f J( , 



i f i f i 

dx dy dz 



u- 2 + A} 2 [u- 2 + Bf 



x < yz — x(3x + y + z) — 



8x 2 \y + z\ 



2x 2 (2yz - 1) 



+3 + 3.x 4 - 2x(l + x z )(y + z) + 



A + B-4 

8(1 - a; 2 ) 2 ] r]9{W) 



A + B 



W 



(29) 



(30) 



ki(k f ) 



(2d) 



+x{y + z){2 - 2a; 2 - yz) + \{Zy 2 z 



Af 



2,2 



y 



, ^3x 4 + x 2 (2 + yz) 



z 



1) 



8x 2 \y + z\ 

2 \-« ~ * ~ -/ ■ A + B _ A 
x(y + z)(3 + 4x 2 yz + 5a; 4 ) - 4 - 3a; 6 - 2a; 4 (l + y 2 + 4yz + z 2 ) 

8(1 -a; 2 ) 2 ] r)9(W) 



+x z (5 - 2y z - 4yz - 2z 2 ) - 



A + B-4 



W 



(31) 



with abbreviations: A = 1 + a; 2 — 2a;y, B = 1 + x 2 — 2xz, W = (x — y) 2 (x — z) 2 — (1 — y 2 )(l — z 2 ) and 
rj = sign[(a; — y)(x — z)]. In fact, only the Pauli-blocked planar box diagram (d) gives rise to a center-of-mass 
tensor interaction from second-order pseudoscalar exchange. 
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Of the three diagrams (e)-(g) in Fig.[5]that encode the effects of medium polarization, (e) and (g) contribute 
to the quasiparticle interaction in the crossed channel and (f) contributes in the direct channel, regardless of 
the form of two-body interaction. In the case of pseudoscalar exchange only (e) and (g) are nonvanishing for 
noncentral components. Diagram (e), representing the coupling of the boson to nucleon-hole states, yields for 
the exchange tensor interaction: 



h (k f ) 



(2e) 



8g 4 M n 



dx 



(1 + Ax 2 ) 4 



2ux 



x 2 ) In ■ 



u — x 



(32) 



MM 



(2e) 



dx 



x 4 (u 2 - 2x 2 ) 



2ux + (u 2 - x 2 ) In 



(1+4x2)4 

The density-dependent vertex correction (g) can be split into a factorizable part: 



u + x 



u — X 



h {k f ) 



(29) 



9 4 M n 



2Air 2 m 3 u 3 



Au 2 



1 + Au 2 



ln(l + 4u 2 



1 + 2u 2 
Au 2 



ln(l + Au 2 ) - 1 



(33) 



(34) 



hi(k f ) 



(2 9 ) _ 



87T 2 m 3 M 5 

and a non- factorizable part: 

MM (V) = - 



1 + 2u 2 
Au 2 



ln(l + Au 2 



dx 



(1 + u 2 )ln(l +Au 2 ) - 3m 2 - 



1 +Au 2 



ln(l +4a; 2 ) - 



4a; 2 



1 + Ax 2 



In 



1 + 4u 2 - Ax 2 
(uVl + Au 2 - Ax 2 + x) 2 ' 



(1 + Au 2 - Ax 2 ) 3 / 2 (1 + Au 2 )(u 2 - x 2 ) 



2ux{l + Au 2 )- 1 



(35) 



(36) 



MM 



(2 3 ' 



9 4 M n 

87T 2 TO 3 M 4 



dx 



ln(l + Ax 2 



Ax 2 



1 + Ax 2 



Aux(l + 2u 2 



(l + 4?i 2 )(l+4u 2 -4a; 2 ) 



-In- 



x t 1 + {u 2 - x 2 )(6 + Au 2 ) (uv 7 ! + 4u 2 - 4a; 2 + x) 2 



(1 + Au 2 - 4a; 2 ) 3 / 2 



(1 + Au 2 )(u 2 - a; 2 ) 



(37) 



We note that diagrams (e)-(g), representing the particle-hole contribution, do not give rise to a cross- vector 
interaction in agreement with the general argument presented in Section ITT] We now evaluate the expressions 
given in eqs. (|2"rj|) -([3"T f choosing g = 5, m = 300 MeV and kf = 1.7 fm. In Table Q] we compare these semi- 
analytical results to those obtained from a numerical evaluation of the second-order contributions as given in 
Section [TT| through a partial-wave decomposition. The agreement is generally on the order of 2% or better. Part 
of this discrepancy (in particular for the Fermi liquid parameter hi) arises from the numerical uncertainty in 
the triple integrations of eqs. (|2"5|) - (|3"TT) . 



B. Spin-orbit interaction to second order 



We consider as well the case of a pure spin-orbit interaction of the form 

2o 2 

V »° = < 2 * fa + 30 ■ (<? X P) » (38) 

[nig + q ) 

where q is the momentum transfer, p is half the incoming relative momentum and m s is the mass of the exchanged 
boson to be fixed later. We consider only the isotropic (L = 0) contributions to the noncentral interactions. 
The first-order contribution from V so to the quasiparticle interaction vanishes trivially since q x p = 0. The 
direct term (a) of the planar box diagram gives the contribution 

with u = kf/m s , while the crossed term (b) contribution reads 

, „ ^(2b) gt M n f 1 , l + 2u 2 1 f, l+4i 2 + 8a; 4 ; „ .1 

h (k f Y 2b > = „f J — In — r+ / da;—- — — - arctan 2a; - arctan x) } . 40) 

v 1! 247rmf\2it 2 1 + u 2 1 + 2u 2 J u 2 (l + 2a;2 3 ^ '] v ; 
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Spin-orbit interaction (kf = 1.7 &n i ) 






[nil J 


2nd(pp+hh) 


0.761 


—0.311 


2nd(ph) 


-0.758 


-0.444 


Total 


0.003 


-0.756 


Analytical (pp+hh) 


0.772 


-0.312 


Analytical (ph) 


-0.761 


-0.448 


Total 


0.011 


-0.760 



TABLE II. The L — noncentral Fermi liquid parameters for the spin-orbit interaction at second-order. The numerical 
results for ho and fco based on the partial-wave decomposition are compared with values from the semi-analytical formulas 
in eqs. p))) (foOTl. 



For the Pauli-blocked planar box diagram, we find that the direct term (c) gives rise to both exchange tensor 
as well as center-of-mass tensor contributions of the form 



h (k f ) 



-^i f l dx fdy /V 



fa 2 k}J J-xJ-x [u- 2 + A] 
Ax 2 \y + z\ 



x(x + y) 



A + B -4 



+ (i-^ 2 ) 2 ( A + B )v <W) 



(A + B- 4)VW 



(41) 



ko(k f ) 



™ = (\ X Cdy \\z 



6n 2 k 3 f J 



A] 4 



x(3x + y) 



3 + 3a; 4 + 2x 2 (2yz - 1) - 2x(l + x 2 ){y + z) + 



Sx 2 \y + z\ 
~ A+B -A 

(1-x 2 ) 2 



A + B- A 



W J 



(42) 



with abbreviations: A = 1 + x 2 — 2xy, B = 1 + x 2 — 2xz, W = (x — y) 2 (x — z) 2 — (1 — y 2 )(l — z 2 ) and 
rj = sign[(a; — y)(x — z)]. Likewise, the crossed term (d) of the Pauli-blocked planar box diagram yields 



h (k f ) 



(2d) _ gsj^n 

6lT 2 k 3 J Q 



dx I dy I dz 



x < yz — x(x + y + z) 



x [u- 2 + A] 2 [u- 2 + B} 2 

Ax 2 \y + z\ +{1 _ x 2 ) 2 (A + B) V 9(W) 



A + B- A 



{A + B -A)VW 



(43) 



k ^ ){2d) - 1% l! dx f ldy f dz 



&x 2 \y + z\ 
A + B -A 



_, { u- 2 + A] 2 [u- 2 + B} 2 

2/o„.~ 1 \ i o i o™4 o^/i i J2 



2x A (2yz - 1) + 3 + 3a; 4 - 2x(l + x 2 )(y + z) + 



yz — x(3x + y + z) 
1-x 2 ) 2 ] 7]0(W) 



A + B- A 



W 



(44) 



The coupling to intermediate particle-hole states (e) in the crossed channel generates both exchange and 
center-of-mass tensor interactions: 



h (k f ) 



(2e) 



9 A s M n 



dx ■ 



(3iru) 2 m 3 s J (1 + 4a; 2 ) 



ux(15x 2 - 17u 2 ) - —{u 2 - x 2 ) 2 In 

2 u — x 



fc (M (2e) 



3ir 2 u 2 m 3 J X (1 + 4a: 2 ) 4 



2 2 \ 1 ^ I X 



2ux + (u — x ) In 



and for the Pauli-blocked crossed box diagram (f), the direct term produces the contributions 



h (k f ) 



(2/) _ *9 4 s M n r 



dx ■ 



(3iru) 2 m 3 Jo (1 + 4a: 2 ) 4 



u(u 2 - x 2 ) + 2u 3 ln- 



Au 2 



„ o i u - x , 3l M + a; 
-bux in — - — - 4a; In ■ 



4a; 2 



(45) 
(46) 



(47) 



fco(fc/) (2/) = 



8g 4 s M n 



dx ■ 



37r 2 M 2 mf Jo (1 + 4a; 2 ) 4 
-2x(u 2 + x 2 ) In 



u(u 2 — a: 2 ) — 4ux 2 In ■ 



4x 2 



u — x 



(48) 
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FIG. 3. Diagrammatic contributions to the quasiparticle interaction in neutron matter generated from the two-pion 
exchange three-neutron interaction. The short double-line symbolizes summation over the filled Fermi sea of neutrons. 
Reflected diagrams of (2) and (3) are not shown. 



Finally, the crossed term of the vertex correction diagram (g) reads 

64 5 fM„ 



dx / dy 



x 3 y 2 



3n 2 k 3 
+ \{u 2 -x 2 ) + 



(1+ 4a; 2 ) 2 (l + % 2 ) 2 



; Re In 



x + \J x 2 



u + y 



: Hey/ x 2 + y 2 — u 2 



xy 



y 



dx I dy 



x 3 y 2 

(1+ 42-2)2(1 + 4y2)2 

2:j; 



:Re\/a 



y 



xy 



(49) 



(50) 



where Re stands for real part. Again, we find that the spin-orbit interaction iterated to second order does not 
give rise to a cross-vector interaction. We choose the parameters g s — 10, m s = 700MeV and kt — 1.7 fin -1 . 
In Table [TT] we compare the sum of the particle-particle and hole- hole diagrams to the sum of diagrams (a)-(d) 
in Fig. [2j Likewise we compare the particle-hole contribution to the sum of diagrams (e)-(g) in Fig. [2j The 
results of both calculations are again in very good numerical agreement with one another. 



C. Spin-nonconserving cross- vector interaction 



As mentioned in subsection III C[ the origin of the spin-nonconserving quasiparticle interaction (proportional 
to the cross-vector operator (cti x a 2) ■ (q x P)) is the interference of the spin-orbit component of the two-body 
potential with its other (central, spin-spin and tensor) components. In order to exemplify this mechanism 
through a solvable model, we consider a simple contact interaction with couplings in all relevant channels: 

Kt = T c + T s Bi ■ a 2 + T t ffi ■ q a 2 ■ q + iT so {a\ + <? 2 ) • (<f x P) ■ (51) 

At second order the interference terms arising from the particle-hole diagrams (e)-(g) in Fig. [2] can be worked 
out analytically. One finds for the first three Landau parameters of the cross-vector interaction 

M n k 3 f 3M n k 3 f llM n k 3 

/ = ^ J -(r c -3r s )r so , h = __J-(r c -3r s )r so , z 2 = - / (3r 8 -r c )r ao , (52) 

57T 707T 847T 

8M n k 5 f 4:M n k 5 f 50M n k 5 f 

io = — ^ — L t r so , h = — — — r t r so , i 2 = — — — L t r so . (53) 

Due to their different structure we have listed separately the interference terms with central and spin-spin inter- 
actions and the "tensor- type" interaction T t c?i • q a 2 • q- Our second-order calculation based on a decomposition 
of the two-body potential into partial-wave matrix elements reproduced these analytical results with good nu- 
merical accuracy. In the absence of the tensor term, the condition T c = 3T S (giving II = 0) is equivalent to a 
vanishing interaction in the spin-singlet state. 



IV. LANDAU PARAMETERS FROM CHIRAL 3N INTERACTION 



In this section we consider the N 2 LO chiral three-nucleon force in neutron matter and derive expressions 
for all L = 0, 1 Landau parameters arising from the leading-order (one-loop) contribution to the quasiparticle 
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interaction. In most cases it is possible to obtain analytical expressions for arbitrary values of L, but for brevity 
we show only the formulas for the isotropic (L — 0) and p-wave (L — 1) parameters. As mentioned in subsection 
III Dt we keep only the two-pion exchange three- neutron force proportional to c\ and c 3 . There are three one- loop 
diagrams contributing to the effective interaction, shown in Fig. [3J which we label as V^^f '* for i = 1,2, 3. Both 
the pion self-energy correction V^Z- and the vertex correction V^^' 2 produce an effective interaction similar 
to that of one-pion exchange. Both terms vanish in the direct channel but give contributions to /, g and h in 
the exchange channel. For the pion self-energy correction we find 

F (k f )^ = (3 - ft ■ ft + 2S n {q))^^ ^-^f - c 3 m 3 + (c 3 - c x )| ln(l + 4m 2 )} , (54) 



F,(k f ) 



(mod,l) 



(3-^-^ + 25x2(9)) 
+ 2(c 3 - ci)u 



g\m\ f (2ci - c 3 )m 



48^ 2 / 4 
3c 3 - 4ci 



2 » 



1 + 4m 2 
ln(l +4m 2 ) 



3ci — 5c 3 )m 



(55) 



with u = kf/m^. As seen from the above formulas, the expressions for the Fermi liquid parameters of the 
central, spin-spin and tensor quasiparticle interaction are identical up to integer factors characteristic of a one- 
pion exchange nucleon-nucleon potential. The contribution from the crossed term of the pion exchange vertex 



correction V, 



mcd,2 



N N 



reads 



Mkf) (mcd ' 2) = (3 - ft ■ ft + 2Si 2 («)) 



(247r)2/^l4« 5 



4m 2 - ln(l + Au 1 ) 
18m 4 + 16m 2 + 3 



c 3 



tr 



m 4 + Au 2 - (1 + Au 2 ) ln(l + 4m 2 ) 
' Au 2 - ln(l + 4m 2 ) ) arctan2M 



+ - 



32m 7 
12u 4 - 16m 6 



ln^l +4m 2 ) + 
30m 2 - 9 



40m 3 „„ 2 
22m + - 

3 M 



3 
2^3 



12m 5 



ln(l + 4m 2 ) 



(56) 



F,(k f ) 



(mod, 2) 



(3 _,V* + 25 12H ))^{^ 



4m 2 - (l + 2M 2 )ln(l+4M 2 ) 



8m 4 + 4m 2 - (1 + 4m 2 ) ln(l + 4m 2 ) 
96m 6 + 80m 4 + 22m 2 + 3 



8 (4m 2 - (1 + 2m 2 ) ln(l + 4m 2 )) arctan 2m 



16m 5 



In (1 +4m 2 ) + 



c 3 

3m 4 , 
56m 6 - 32m 8 



60m 4 



48m 2 - 9 



6m 3 



ln(l + 4m 2 



28m 3 



3 

10m + - 

M 



(57) 



Lastly, we compute the L = 0, 1 Landau parameters associated with the Pauli-blocked two-pion exchange 



diagram V^r d ' 3 , which has a more complicated spin structure than V^ 1 ^ 1 ' 1 and V } 



d.l 



od.2 



NN 



Both the direct and 



exchange terms contribute, and the central components of the quasiparticle interaction read 



F (k f ) 



(mcd,3) 



9 2 A m l 
16^ 2 / 4 



(C 3 - CljM 



8c 3 m 3 3c 3 - 4ci 



-(12ci — 10c 3 ) arctan 2m + (1 + ft ■ ft) 



u 

1 

dx 



ln(l + 4m 2 ) 



C 3 , 



2ci Z z + ^-(X 2 + 2Y 2 ) 



^ l(fc/) (mcd, 3) = (1 + 3i . *? a )^| [dx{*cv& + Ml) + c 3 (X 2 a + 2X 2 + \x< 



(58) 



(59) 



The auxiliary functions X, Y,Z; X a , X} }1 X c ; and Z a ,Zt, are defined in Section 2.2 of ref. [2j. Concerning the 
noncentral quasiparticle interaction, V^^' 3 produces no tensor forces in pure neutron matter. However, one 
expects V^^f' 3 to generate an exchange tensor force in symmetric nuclear matter [32j since in this case the three- 
nucleon force proportional to C4 does not vanish. As a special feature, the crossed term of the Pauli-blocked 2ir 
exchange diagram gives rise to a cross-vector interaction. With its usual representation given by 



Across = (ft x ft) • (q x P) ^2 h{k f ) P L {Pi ■ P2) , 



(60) 



L=0 
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a complete analytical solution for the Landau parameters Zz(Jfc/) could not be obtained. When choosing the 
alternative form of the cross vector interaction 



•Across = — . ^ |2 2^ lL ( k f> P LKPl -P2), 

L=0 



Wl +P2| 



(61) 



all occuring integrals can be solved analytically in the present case. The results for the Landau parameters read 



UkA( med >V - 9Am ^ — 



+C 3 



4?r - 8w h 



16u 4 - (1 + 4u 2 )ln 2 (l + 4u 2 ) 
1 1 + 4m 2 



u 2u 3 



ln(l + 4u 2 



1 1 



16u 5 



ln 2 (l + 4m 2 ) 



(62) 



Wf) - (87r)2/ 4^ c i 



8u- 



32u 6 + 24u 4 + 8u 2 + 1 



2 8u 4 + 6u 2 + 1 
— + 5 

ln 2 (l + 4?i 2 ) 



ln(l + 4u 2 



+c 3 



20u 3 2 3 3 16u 6 + 28u 4 + 18m 2 + 3 
16u j ? H 5 

3 u u A 4w tl 8u' 



x ln(l + 4u 2 ) + 7-r-?r(64u 8 - 20u 4 - 8u 2 - 1) ln 2 (l + 4u 2 ) 
64u J 



(63) 



The numerical results assuming the standard form, eq. (|60[) , of the cross vector interaction will be shown in the 
following section. 



V. RESULTS 



A. Neutron matter equation of state 



The Fermi liquid parameters in neutron matter, unlike those in symmetric nuclear matter close to the sat- 
uration density, are largely unconstrained by empirical data. Given that our perturbative treatment of the 
quasiparticle interaction includes only the leading-order medium corrections from two- and three-body forces, 
it is useful to compare the associated zero temperature equation of state of neutron matter (at the same order 
in perturbation theory) with those obtained using nonpertubative methods. As a benchmark we consider vari- 
ational calculations [34j of neutron matter employing the high-precision Argonne v\% two-nucleon potential [35} 
together with the Urbana UIX three-body potential [36( , which provide a realistic description of light nuclei 
and nuclear matter. 

In the partial-wave representation of the two-body interaction the first-order contribution to the energy per 
particle E — E/A reads 



^ ) (fc/) = ^-43E( 2J+1 ) dpp 2 (k f -p) 2 (2k f +p)(plSJ\V\plSJ) 



(64) 



V isj Jo 
and the more intricate second-order contribution takes the form 

'0 



(47r) 4 fc3 



E 



Sm B m'JJ'M 



2 /-V fc /-P"/4 

dp' p' / dpp 2 







dq q 2 N (limhm' Ismhm') 



x 1 



p 2 — q 2 



d cos QpP™ (cos 6p)P™ (cos 6p) / dcos6 q P% (cos6 q )P% (cos9 q ), 



(65) 



where x p = min{l, (k 2 — p 2 — p' 2 / '4) / 'pp'} and x q — min{l, (q 2 — k 2 + p' 2 / '4) / qp'} . In the Appendix we provide 
the analytical expressions for the second-order contributions to the energy per particle associated with the two 
model interactions introduced in Section [IIII Finally, the leading-order chiral three-neutron interaction leads to 
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FIG. 4. Energy per particle of neutron matter from chiral and low-momentum two- and three-body interactions. The 
cutoff scale associated with the bare chiral nuclear potential is A = 2.5 fm -1 , while that of the low-momentum interaction 
is A = 2.1 fm -1 . The curve labeled 'APR' is taken from the variational calculations of Akmal et al. [§3], and the curve 
labeled 'QMC is taken from the quantum Monte Carlo calculations in refs. 39, 4(j. 



Hartrce and Fock contributions to the energy per particle of neutron matter of the combined form 

' 3 2c 3 



„ . (6ci - 5c 3 )^- arctan2u 



"3„.6 



u b + (c 3 - Cl )u 4 



+ (3 Cl -2c 3 )— + 



i2-¥ + T (3c3 " 4ci) 



ln(l + 4u 2 ) 



32u 3 J 



dx 



6 Cl H 2 + c 3 (G| + 2G|) 



(66) 



with u = kf/m^. The auxiliary functions H, G5 and Gt are defined in eqs. (24-26) of ref. [37[. These 
interaction contributions are added to the relatististically improved kinetic energy per particle E^{ n {kf) = 
3k 2 f /10M n -3kj/56Ml 

In Fig. |4] we plot the resulting equation of state of neutron matter for densities up to p ~ 1.5y0o- For both 
chiral and low-momentum interactions one finds good agreement with the results of the variational calculation of 
ref. [jy]], labeled 'APR' in the figure. The difference between the two perturbative calculations arises primarily 
from the different values of low-energy constants 01,3. Note that the recent neutron matter calculation [38[ 
including the subleading chiral three- and four-neutron interactions gives very similar results. For comparison, 
we have included in Fig. 0] the result for the neutron matter equation of state obtained in a recent quantum 
Monte Carlo (QMC) calculation [39| employing phenomenological density-dependent two-body potentials. The 
data points for the lowest two densities are taken from a different quantum Monte Carlo calculation in ref. [4fj| . 



B. Neutron matter quasiparticle interaction 



In this section we present and discuss the calculations of the L — 0, 1, 2 Landau parameters for the quasiparticle 
interaction in neutron matter from chiral and low-momentum two- and three-nucleon forces. We begin by 
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-0.112 


-1.188 


-0.679 
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1.053 


0.613 


0.337 


0.197 


1.212 


0.654 


0.346 


0.195 


h [fm*] 


0.270 


0.060 


-0.040 


-0.080 


0.239 


0.102 


-0.051 


-0.079 



TABLE III. The L = 0,1, 2,3 Fermi liquid parameters of the bare N 3 LO chiral NN potential of ref. 19] as well as the 
low-momentum NN potential Viow-k at a resolution scale of 2.1 fm -1 at first order in many-body perturbation theory. 

considering the leading-order (free-space) contribution from both the bare chiral N 3 LO NN interaction as 
well as the low- momentum NN potential Mow-k [HI obtained by integrating out momenta beyond the cutoff 
scale of 2.1 fm -1 . We use the general formula in eq. (j9|) for computing the first-order perturbative contribution 
to the quasiparticle interaction as well as the projection formulas in eq. (|7|) for extracting the scalar functions 
/, g, and h. In Table Hill we show the results for neutron matter with a Fermi momentum of kf — 1.7 fm 
corresponding to a density of po = 0.166 fm -3 . In both cases the Fermi liquid parameters decrease rapidly in 
magnitude with L for all channels. For larger values of L the difference between bare and evolved two-neutron 
interactions is strongly reduced, and at L = 3 it is almost negligible. Short-distance repulsion in the bare chiral 
NN interaction is integrated out through the renormalization-group evolution, leading to a significant reduction 
in fo- For both potentials, the compressibility K. of neutron matter at po would be unphysically small [l4l ]. The 
isotropic component of the spin-spin interaction, go, is enhanced by ~ 15% at the resolution scale of 2.1 fin. , 
while the isotropic component ho of the exchange tensor interaction is reduced by a slightly smaller factor. 
The p-wave component of the spin-independent quasiparticle interaction, fx, increases as the resolution scale 
decreases. The corresponding values of the quasiparticle effective mass are M* /M n = 0.78 and 0.84 for the 
chiral and low-momentum interactions respectively. 

Next we compute the second-order particle-particle, hole-hole and particle-hole diagrams, shown in Fig. [TJ 
with two-neutron interactions. These provide the leading-order Pauli-blocking and polarization effects and 
give rise to components of the quasiparticle interaction depending explicitly on the center-of-mass momentum 
P = pi +P2- In Table |IV] we display the Fermi liquid parameters associated with the second-order contributions 
(j>p, hh and ph) for both the chiral N 3 LO nucleon-nucleon interaction and the low-momentum interaction 
V|ow-k in neutron matter at a density po- For the isotropic components, it is generally true that the particle- 
particle diagram gives contributions that are significantly larger than the hole-hole diagram. However, for 
higher values of L this is not the case, and the hole-hole diagram is in general comparable in magnitude to the 
particle-particle contribution. Coherent effects among the three diagrams arc observed especially for the Landau 
parameters fx, go, ho and &q. Second-order effects thus tend to dramatically increase the quasiparticle effective 
mass M* , decrease the spin susceptibility of neutron matter and reduce the isotropic exchange tensor strength 
ho. Although individually large, the contributions to fo nearly cancel at second order for the bare chiral NN 
interaction. With the low- momentum interaction the sum gives fo — 0.379 fm 2 , which approximately cancels 
the reduction in fo from the renormalization group evolution at first order. In general, the particle-particle 
contributions decrease significantly in magnitude as the resolution scale is lowered. This effect is due primarily 
to the reduction in phase space in the particle-particle channel as the momentum cutoff is lowered. The results 
compiled in Table ITVl were obtained with free particle energies in the denominators of eq. (fT2")l . In all subsequent 
tables, we include as well the one-loop corrections to the single-particle energies in the second-order diagrams 
(for details see ref. [l[). 

Finally, we calculate the contributions to the Fermi liquid parameters from the leading-order chiral three- 
neutron force. We evaluate the analytical formulas in Eqs. (l54l) - ([59|) for the central and exchange tensor 
contributions and perform numerical calculations of the cross- vector Fermi liquid parameters lh{kf) m the 
standard representation, eq. (|60l) . In Fig. [5] we plot the density-dependent Fermi liquid parameters for the 
chiral three-neutron force with low-energy constants cx = —0.81 GeV -1 and C3 = — 3.2GeV _1 . For densities 
larger than p ~ 0.5po the Landau parameters scale approximately linearly with the density. The largest effect 
is a strong additional repulsion in the isotropic spin- independent parameter fo- In fact, at nuclear matter 
saturation density po, the strength of the three-body correction in this channel is larger than the two- neutron 
force contributions at 1st and 2nd order together. The quasiparticle effective mass M* , governed by the 
parameter fx , is reduced by less than 5% at saturation density po with the inclusion of the two-pion exchange 
three-neutron force. Similar observations have been made in ref. @ for the case of symmetric nuclear matter 
where additional three-nucleon forces proportional to the low-energy constants C4, cd and ce are present. The 
qualitative similarities between the quasiparticle interaction in neutron matter and symmetric nuclear matter 
from three-body forces are due to the dominant role played by contributions proportional to the low-energy 
constant C3 [2j. Neglecting effects from noncentral Fermi liquid parameters, the chiral three- neutron force at 
leading order tends to enhance the spin susceptibility x of neutron matter (see eq. ©) by lowering go by 
about 25% at po from the value obtained with two-nucleon forces only. The three-body correction to the 
exchange tensor interaction is smaller in magnitude than the central force contributions. For the exchange 
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-0.090 


-0.066 
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-0.044 


Total 


-0.212 
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-0.143 


-0.192 


-0.012 


-0.080 
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PP 


-0.085 


0.064 


0.014 
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fc [fm 2 ] 




hh 


-0.036 


0.052 


-0.017 
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fc [fm 2 ] 




ph 
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-0.017 


0.075 
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Total 


-0.178 
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0.072 


-0.119 


0.056 


0.043 
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ph 


0.135 


-0.031 


-0.279 


-0.062 
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-0.161 



TABLE IV. Second-order contributions to the L — 0, 1, 2 Fermi liquid parameters in neutron matter characterized by the 
Fermi momentum kf = 1.7 fin. . We have separately listed the particle-particle (pp), hole-hole (hh) and particle-hole 
(ph) contributions for both the bare N 3 LO chiral NN interaction as well as the low-momentum NN potential VL>w-k with 
Ai ow _ k = 2.1 fm" 1 . 
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tensor interaction the pion-self energy correction V^ff' 1 and vertex correction ' 2 enter. As seen already in 
ref. [2|, |32j both of these terms have the structure of one-pion exchange and are individually large with opposite 
sign. The cross- vector interaction from chiral three-nucleon forces is relatively small compared with the central 
quasiparticle interactions. However, at saturation density p Q the magnitudes of the cross-vector Fermi liquid 
parameters /o,i,2 are comparable in magnitude to those from two-body potential at second order. 

In Fig. [5] we combine the results for the first- and second-order two-body contributions with the leading- 
order three-body corrections at a resolution scale of A = 2.5 fm -1 . All Fermi liquid parameters up to L = 2 
are plotted as a function of the neutron density p. The same quantities are compiled in Table [V] at a Fermi 
momentum of kf — 1.7 fm -1 together with the corresponding values from the low- momentum interaction VJ ow _k 
at the scale A = 2.1 fm -1 . From Fig. [5] one observes that all of the Fermi liquid parameters considered here 
for the spin- independent part of the quasiparticle interaction vary strongly with the density up to p ~ pq. In 
particular the parameter /o has a nearly linear depedence on the density, and at kf = 1.7 fm -1 gives rise to 
a compressibility JC — 560 MeV of neutron matter that grows strongly with increasing density. The Landau 
parameter f\ decreases rapidly at small densities, but flattens out close to nuclear matter saturation density, 
where it nearly vanishes thereby yielding an effective mass M*/M n ~ 1 according to eq. (J5J. The L = 2 
component of the spin-independent quasiparticle interaction exhibits a nearly linear decrease with the neutron 
density owing to coherent effects from second-order two-nuclcon forces and three-neutron forces at leading-order. 
The isotropic component of the spin-spin quasiparticle interaction, go, decreases with the neutron density, but 
we find no evidence for a phase transition to a ferromagnetic state (Go = N g < —1) close to nuclear matter 
saturation density. For densities larger than p , the density-dependence of go is governed mainly by the three- 
neutron force contributions, which decrease go and consequently increase the neutron matter spin susceptibility 
X- The isotropic component ho of the exchange tensor interaction is nearly independent of the density, since 
the medium corrections from two- and three-body forces approximately cancel. The Fermi liquid parameters of 
the novel center-of-mass tensor and cross- vector interactions are relatively small in magnitude, but ko, lo and 
I2 have a strong density dependence. From Table [V] we see that there remains a moderate dependence on the 
resolution scale and choice of low-energy constants. These variations lead to differences on the order of 10-30% 
for most of the L = 0, 1, 2 Fermi liquid parameters. 
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TABLE V. Fermi liquid parameters for the quasiparticle interaction in neutron matter at a density corresponding to 
a Fermi momentum of kf = 1.7fm _1 . The low-energy constants of the N 2 LO chiral three-nucleon force are chosen 
to be Ci = —0.81 GeV -1 and C3 = —3.2 GeV -1 when employed together with the bare chiral NN interaction and 
ci = —0.76 GeV -1 and C3 = —4.78 GeV -1 with the low-momentum interaction Viow-k- 



VI. CONCLUSIONS AND OUTLOOK 



In the present work we have computed the L = 0, 1, 2 Fermi liquid parameters for the quasiparticle interaction 
in neutron matter employing realistic two- and three-nucleon interactions derived within chiral effective field 
theory. In addition to the free-space contribution from the two-body interaction, we have calculated without any 
simplifying approximations the second-order two-body corrrection as well as the leading three-body correction. 
A general method for extracting all components of the quasiparticle interaction, both central and noncentral 
parts, for two-body interactions given in a partial-wave representation has been developed and tested with 
simple model interactions that can be treated semi-analytically. Employing realistic two- and three-neutron 
forces, we find that medium-dependent loop corrections play an important role in increasing the compressibility 
of neutron matter from an unphysically small value to about K, = 600 MeV at nuclear matter saturation density 
Pq. Second-order effects from two-body forces strongly enhance the quasiparticle effective mass M*, while 
three-neutron forces play only a minor role for this quantity. The first- and second-order contributions to 
go from two-body forces are positive for the densities considered in the present work, though they decrease 
in magnitude for increasing density. When combined with the medium-dependent loop corrections from the 
leading-order chiral three-neutron forces, the Landau parameter go decreases with density, which leads to an 
increasing spin susceptibility x of neutron matter. The noncentral components of the quasiparticle interaction 
Hl, fci and II have been computed as a function of the neutron density p. The extent to which they affect the 
spin susceptibility of neutron matter as well as the response functions for electroweak probes will be studied in 
future work. 



VII. APPENDIX: ENERGY PER PARTICLE FROM MODEL INTERACTIONS AT SECOND 

ORDER 

Organizing the second-order calculation of the energy density pE in the number of medium insertions, see 
eq. (|22[) . one must evaluate (two-ring) Hartree and (one-ring) Fock diagrams each with two or three medium 
insertions. We give first the pertinent analytical expressions for the pseudoscalar interaction in eq. (|23p at 
second order. 

Hartree diagram with two medium insertions: 

E{k f )^ = - 15k + 32arctan2u - -^(3 + 20u 2 ) ln(l + Au 2 )\ , (67) 

(87r) cl L 2u Su a J 

Fock diagram with two medium insertions: 

E(k f ) (F2 ' ) = (4 ~% j 2 " 2 ~ 7u3 arctanu + 6V2u 3 arctan(V2u) 

+ (l + ^) ln(l + u 2 ) - Z(l + 2u 2 ) ln(l + 2u 2 ) 

. f\ / n2/„ x 1 + 4a; 2 + 8x 4 . 

+4 / dxxiu — x) (2u + x) — ; tttt; — (arctanx — arctan2x) 

Jo (l + 2x 2 ) 3 

Hartree diagram with three medium insertions: 
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(69) 
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Fock diagram with three medium insertions: 
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For the spin-orbit interaction, eq. ([55]). at second order the analogous expressions read 
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(70) 



(71) 



(72) 



(73) 



(74) 



where u = kf /m s . An interesting feature of the second-order spin-orbit interaction is that for large mass m s 
the Hartree and Fock contributions become equal with the same sign. It is worth noting that the agreement 
between these semi-analytic results and those based on the partial wave decomposition eq. (|65[) agree on the 
per mille level. 
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